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Abstract 

The study of dilepton and direct photon emission was one of the main topics of the experimental 
program at the SPS devoted to the search of signals for QGP formation. Three generations of 
experiments, Helios-3, NA38/NA50, CERES and NA60 measured <? + e~ or ii^yT production in 
various colliding systems and at different energies. While lepton pair production in p-A collisions 
was found to be reasonably well described by the expected sources, all experiments observed in 
nuclear collisions an excess of the yield above the extrapolation from p-A. As a result of this 
joint experimental effort we have currently a large amount of information characterizing this 
excess: its mass spectrum over the full range from 0.2 GeV/c 2 up to the //i/f, its transverse 
momentum spectra including their mass dependence, its angular distributions, its dependence on 
collision centrality over the complete range etc. Putting together all this information leads to the 
conclusion that what we observe is the long-sought thermal radiation from the fireball. 



1. Introduction 

The properties of the dilepton and real photon production in relativistic heavy ion collisions 
constitute a significant fraction among the various signatures of the Quark-Gluon Plasma (QGP) 
formation. Being produced at all stages of the collision and practically lacking the final state 
interactions, they convey dual information. From one side, their production rate and kinematics 
are sensitive to the properties of the surrounding matter: dominance of the partonic or hadronic 
constituents, temperature, density, flow etc. In turn, the hot and dense matter is expected to affect 
spectral functions of the dilepton emission. Disentangling these two effects poses both inter- 
petational and experimental problems: the "conventional" sources of the observed effects need 
to be understood, and the small production rates and large backgrounds require large integrated 
luminosity measurements. 

2. Experimental results 

The mentioned difficulties are the reason behind the scarce results on direct photons collected 
so far at the SPS. 15-20% upper limits (95% CL) on the excess over hadronic sources tj and 
j]' decays) were set by WA80 HI and in the central S-Au collisions at 200 A GeV. The most 
significant result [3| is the observation of up to 20 ± 7% excess for pj > 1.5 GeV/c in Pb-Pb 
collision at 158 A GeV. It is well described by theoretical models involving pQCD photons and 
thermal emission both from hadron gas and QGP [4|, but the interpretation is ambiguous: the 
thermal photon emission rate is determined by the temperature and close to the transition point 
both phases provide similar rates. 
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The situation is much more advanced in the dilepton sector. Conventionally, its mass spec- 
trum is separated into three regions, distingished by their dominant contributions in p-p and p-A 
collisions: (i) Low Mass Region (LMR, M < M^) where the emission is mostly determined by 
the leptonic decays of vector mesons (resonant part) and non-resonant contribution from Dalitz 
decays; (ii) Intermediate Mass Range (IMR, < M < Mj/$) composed of Drell-Yan dileptons 
and pairs coming from uncorrected decays of open charm (mostly D and D mesons); (iii) High 
Mass Range (HMR, M > Mj^) dominated by the heavy quarkonia and Drell-Yan pairs. This 
part of the spectrum was covered by in this volume and will not be considered here. 

2.1. Excess in LMR and IMR 

With arrival of the ion beams on CERN SPS all dilepton experiments observed an enhance- 
ment of the l + r production with respect to the extrapolation from p-A collision. Helios/3 |6| 
and NA38 Q experiments reported the continuum excess in fi + ^r , comparing the S-W and S-U 
collisions at 200 A GeV respectively with p-W interactions at the same energy. 

An important milestone was set by the CERES experiment which measured electron pairs us- 
ing a system of two RICH detectors. It observed a strong (5. ± Q.l(stat) + 2.{syst)) LMR excess 
in e + e~ production in S-Au collision [8| with respect to the "cocktail" of hadronic decays [9|, 
well describing the p-Be and p-Au spectra [9|. The triggered theoretical activity, trying to de- 
scribe these data, was not successful until the appearence of calculations involving in-medium 
7T7T annihilation via a modified p spectral function. The model [10], exploring the idea of the 
Brown-Rho scaling [11], assumed a decrease of the p pole mass close to the chiral symmetry 
restoration point, while lfl2l predicted a brodening of the p, with significant contribution from 
the interactions with baryons. Unfortunately, the low statistics and insufficient (~ 6%) mass res- 
olution did not allow to give preference to one of the models. The excess was later confirmed 
(2.73 + 0.23(stat) ± 0.65(syst) ± 0.82(decays)) in Pb-Au interactions at 158 A GeV. Fig^left) 
shows the e + e~ combined 95/96 spectra together with the calculations ifTOl and lfT2l . 




M ee (GeV) m ee (GeV/c 2 ) m ee (GeV/c 2 ) 



Figure 1: CERES inclusive e + e mass spectra for Pb-Au at 158 A GeV compared to models |10| and |12 |. (Left): 
combined 95/96 data 1 13), (Center) 2000 da\a,(Right) Same data after the subtraction of the hadronic "cocktail" 1141 . 

After the CERES upgrade by the TPC in 1999 (which improved the mass resolution to ~ 4%), 
the only existing SPS dilepton measurement at 40 A GeV was done [15|. Enhancement 5.9 ± 
\.5{stat) + \.2{syst) + \&{decays) was reported in Pb-Au collisions, stronger than at 158 GeV. 
The latter was remeasured in 2000 [14| with the result 2.6 ± 0.3(sfaf) + 0.4(^5?) + 0.8(decays), 
see Fig.^center). Fig.^right) shows the isolated excess of 2000 Pb-Pb data, with the hadronic 
decay contributions subtracted using the rates provided by the statistical model. The possible 
explanation of the significantly larger excess at lower energy is the higher baryonic density due 
to the stronger stopping. 
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In the same period another important result was provided by NA50 ifTBI . Using the muon 
spectrometer inherited from NA38 and equipped with a very selective pp. trigger system, it ob- 
served in Pb-Pb collisions at 158 A GeV an IMR excess with respect to the expected Drell-Yan 
and open charm contributions. These latter were extrapolated from the p-A dimuon spectra at 
450 GeV. Combined with the reanalyzed S-U data of the NA38, the excess showed an ap- 
proximately linear rise with the number of participants, reaching a factor 2 in the central Pb-Pb 
collisions. The excess, which had a kinematics resembling that of the open charm, could be 
accounted for IfTTl either by some mechanism enhancing the latter contribution in the NA50 
acceptance window or by the long sought thermal dimuons 1 18 1. 

The recent results are dominated by the precise and high statistics dimuon measurements 
of the NA60 experiment [19]. Its upgrade of the NA50 setup by a radiation tolerant silicon 
pixel vertex tracker (VT) [20 1 placed in a 2.5 Tesla dipole magnet between the target and the 
hadron absorber provided numerous advantages with respect to its predecessors. Particularly, by 
matching the muons from the spectrometer to tracks reconstructed in the VT, NA60 significantly 
improved the dimuon mass resolution (from ~ 80 MeV/c 2 to ~20 MeV/c 2 at u> mass) and was 
able to measure the offset of the muons wrt the production vertex with a resolution of ~50 pm: 
enough to distinguish between prompt dimuons and those coming from the open charm decays. 

In 2003 NA60 has collected ~230 millions dimuon triggers from In-In collisions at 158 
A GeV . Data were taken with two current settings in the spectrometer magnet: 4 kA for enhanced 
acceptance at low mass and pj and 6.5 kA to enhance the IMR and HMR statistics. 
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Figure 2: NA60 dimuon mass spectra measured in In-In collisions at 158 A GeV: (Left) peripheral |21 1 (Center) inte- 
grated over centrality (excess is shown by triangles) |22 | (Right) Isolated excess in semicentral collisions [ 22'] compared 
to the Brown-Rho and Rapp-Wambach models 1 23 1. The "cocktail" p is not subtracted. 

Figpj'Ze/y shows the signal dimuon mass spectrum for peripheral (dN c h/dr] < 30) In-In colli- 
sions 1211 . It is well described by the "cocktail" of hadronic decays obtained using the GENESIS 
generator [9| improved and adapted for dimuons [24 1 . The high statistics and good mass reso- 
lution allowed direct fit of the rj — > p* p~y and a> — > p + p~n° decays form-factors ll25l with the 
pole approximation F = (1 - M 2 /A 2 Y\ yielding A~ 2 (in GeV' 2 ) 1.95±0.17(stat.)+0.05(syst.) 
and 2.24±0.06(stat.)±0.02(syst.) respectively. The values agree with previous measurements 
by the Lepton-G experiment 11261 . improving their errors and confirming the strong enhance- 
ment of the u> form-factor with respect to the VMD expectation of A~ 2 = 1.68 GeV~ 2 l27l . 
These measurements, together with the improved value of the u> — > p + p~jr branching ratio 
([1.73+0. 25(stat.)±0.14(syst.)]-10~ 4 ), significantly decrease the uncertainty in the hadronic de- 
cay "cocktail" below the a> in the analysis of more central collisions. Fig^center) shows the 
LMR mass spectrum integrated over centrality [22|. Thanks to the good mass resolution the 
narrow peaks of 0, u> and rj as well as the rj, r\' and ai Dalitz decays can be locally subtracted 
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uncovering a significant excess (shown by triangles) centered around the p (its "cocktail" contri- 
bution is not subtracted). Figj2j>/g/if) compares the isolated excess for the semi-central collisions 
(dN c h/dT] = 140) (22) with the preditictions of Brown-Rho [10] and Rapp-Wambach HH mod- 
els ll23ll . The complete disagreement with the dropping mass and nice agreement with the broad- 
ening scenario is obvious. A similar conclusion fl4l . although with less statistical significance, 
is obtained by CERES from the Pb-Au (2000) data (Fig.^right)). 
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Figure 3: Evolution of the LMR excess with the multiplicity. (Left) NA60, In-ln at 158 A GeVEfl. (Right) CERES, 
combined 95/96 Pb-Au at 158 GeV data (B). See text for details. 

Fig.[3]7e/f) shows the centrality dependence of the LMR excess studied by NA60. The total 
"excess" (circles) is separated into continuum (filled triangles) and peak (hollow triangles) parts 
and related to expected "cocktail" p contribution (assuming p/oj — 1 on the cross section level). 
While the peak contribution agrees with the p produced at freeze-out, the continuum part shows 
a monotonic rise and broadening with centrality (21]. Such a faster than linear rise of the excess 
with multiplicity is compatible with emission from the annihilation process. It can be conjec- 
tured that the magnitude of the continuum excess directly measures the lifetime of the fireball in 
number of p generations: so called p-clock [28 1. The centrality dependence of the excess seen 
by CERES [13] in Pb-Au collisions, shown in Fig. [SJng/tfJ as a ratio of the excess to the total 
"cocktail" expectations, is in a good agreement with NA60 results. 

The analysis of the NA60 IMR mass spectrum |29| in the In-ln data is done by fitting the 
signal in the 1.16 < M < 2.56 GeV/c 2 range with a superposition of the Drell-Yan and open 
charm contributions obtained from the Pythia 6.325 generator. The fits are done in terms of 
multiplicative factors for the reference cross sections. The latter are defined in the following 
way: for the Drell-Yan, it reproduces the cross sections measured above the if/' by NA3 [ 30 1 and 
NA50 [ 3 1 1, while for the open charm the result of a similar fit to NA50 p-A data at 450 GeV ifTBI 
rescaled to 158 GeV by Pythia is used (<x cf = %.6p.b). Due to the insufficient high-mass Drell- 
Yan statistics, the integrated effective luminosity is extracted from the number of Jj\p events 
and its cross section (corrected for the nuclear and anomalous suppression effects). Fig. ^left) 
shows such a fit to the centrality integrated dimuon mass spectra with low current setting in the 
spectrometer. At this level the results are fully compatible with observations of NA50 in Pb-Pb 
data ifTBTl : while the Drell-Yan contribution agrees well with the expectations (and the data above 
if/'), a strong excess with a mass shape resembling the open charm contribution is found. The 
global fit to both data sets leads to the enhancement factors 1.26 ± 0.09 for Drell-Yan and 2.61 ± 
0.20 for open charm. To clarify the origin of this excess NA60 uses its excellent muon offset 
resolution to separate (statistically) the open charm (off-vertex decays) and prompt contributions. 
Fig. fright) shows the dimuon "offset'^ distribution for the mass range indicated in the left 



'Defined as A,, = (dV 'd T )'' 2 for the single muons and A w 

4 



[(A„i + A„2)/2] 1/2 for the dimuons, with d and V 



panel, fitted to the expected open charm and prompt (assumed to be Drell-Yan) contributions. 
The global fit to low and high current data sets leads to enhancement factor of 2.29 + 0.08 and 
1.16 ±0.16 for the prompt and open charm samples. Hence the excess should be attributed to the 
prompt emission, while the open charm contribution is compatible with the extrapolation from 
the NA50 pA data. 




Figure 4: IMR signal dimuon distributions in In-In collisions at 158 GeV (29 1 : (Left) mass spectra fitted in IMR by the 
superposition of Drell-Yan and open charm contributions; (Right) same for weighted offset spectra. 

The statistics of the NA60 IMR spectra is not high enough to extract the open charm con- 
tribution differentially (in centrality, M, pj etc.). For this reason, in the differential analysis 
the kinematic distributions for open charm and Drell-Yan are taken from Pythia spectra, while 
for the dependence on centrality both are assumed to scale with the number of binary collisions 
extracted from the measured number of Jj^i events in a given dN c i,/drj bin (corrected for the 
suppression). The IMR excess is defined as the difference between the measured signal and the 
sum of these two contributions. 

The dependence of the excess on the number of the collision participants shows a faster than 
linear scaling with the number of binary collisions (with excess/DY reaching a factor 2.3 ± 0.7 
for the most central with N par , > 200), but is slower than quadratic increase with the squared 
number of participants! 29 1. As in the case of the LMR excess such a behaviour is compatible 
with the emission from the annihilation in thermalized medium. Fig. [5]7e/f) summarizes the 
mass spectrum of the excess seen in the In-In collisions ||29l , corrected for the acceptance and 
reconstruction efficiency and normalized to the per charged particle yield. 

NA60 performed also an extensive study of the pp. mj spectra l32l . summarized on Fig.^right). 
Both hadrons and excess mj spectra are well described by the dN/dnij oc exp(-mj /J ' e jf) form, 
except for some puzzling softening of the excess at masses < 1.2 GeV/c 2 in the range of 
mj - M < 0.2 GeV . All hadrons - rj, a>, p (defined as the peak on top of the LMR excess) 
and (p show the rise of T e g with mass characteristic for the "blue-shift" due to radial flow. The 
deviations T p > Tq, from nearly linear scaling with mass (T e g m Tq + M < f3 > 2 ) are compat- 
ible with different freeze-out times of the corresponding hadrons due to their different coupling 
to the expanding medium. This is supported by a "blast-wave" analysis 0331 : the (f> decouples 
first, when the flow is not yet fully developed, while the p, whith its strong coupling to pions, 
freezes out last and profits from the full flow (with T e g reaching ~ 300 MeV in the most central 
collisions). A similar increase of T e ff with mass is observed for the LMR excess. Since, due 
to the "soft point" in the equation-of-state [34| the (eventually produced) partonic phase at SPS 



being the vector and corresponding covariance matrix of the transverse offset of the muon wrt. the primary vertex 
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energies is not expected to develop significant flow, this suggest an emission from the thermal- 
ized hadronic expanding gas. Surprisingly, the rise of the excess T e fi with mass changes to the 
flat behaviour after the sudden drop by nearly 50 MeV at M ~ 1 GeV/c 2 . The most plausible 
explanation of this effect is that the excess at these masses is dominated by thermal dimuons 
from a partonic source lacking significant flow. 
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Figure 5: (Le//)Acceptance corrected mass spectrum of the excess (the "coctail" p is not subtracted) in In-In 

collisions | 29 | (Right) T c ff of excess (the p peak contribution is subtracted) vs. dimuon mass f33] |35l . 

Such an interpetation is supported by the theoretical calculations. The models labeled Hees/ 
Rapp IT361 and Renk/Rupppert [37] assume the dominant contribution in the LMR to be n + 7T~ 
annililation via modified p spectral function ([12] and [38 1 respectively), while the IMR excess 
is defined by the combination of 4 n annihilation processes and a significant contribution from the 
annihilation in the partonic phase. They differ in the estimate of the fraction of the latter: ll36l puts 
emphasis on chiral mixing via the na\ — > p + pT process, keeping the partonic contribution within 
20 - 60% (depending on the fireball evolution scenario), while 071 assumes ~ 80% contribution 
from the partonic phase. The model Dusling/Zahed [39] uses hydrodynamic calculation with a 
virial expansion for the rates in the hadronic phase and qq annihilation in the partonic one. The 
latter contributes 60 - 90% to the IMR excess. All models roughly agree with the data. The 
differences at low masses reflect the differences in the tail of the p spectrum ( ||36ll with a strong 
effect from baryons provides the best description). A detailed comparison of these models with 
data in narrow pj bins is contained in [33]. 

The thermal origin of the excess is further supported by the abscence of any polarization in the 
excess dimuons [40 41 j. All coefficients of the dcr/dQ. cc (l+Acos 2 6+psin28cos(p+^sin 2 0cos2(p) 
parameterization are found to be compatible with zero, which is a necessary (though not suffi- 
cient) condition of the emission from an isotropic thermalized source. 



2.2. u> and cf) mesons 

Despite various predictions (|36[ and references therein), no anomaly in the pole positions 
and widths of the u> and (f> mesons was observed so far within the experimental reach of NA60. 
Due to the longer lifetimes, only a small fraction of these mesons produced in the medium 
contributes to the observed dilepton signal. For this reason any in-medium modifications will 
have much weaker effect on the measured u and (f> spectra than for the p. Besides that, the 
possible broadening of the a> would be practically inobservable since it would merge with the p 
peak. Instead one could look for the deficit of low-pr cl> dileptons (with strongest contribution 
from in-medium decays) in the nominal pole position. NA60 has reported the first observation 
of such an effect [32 1 for the to. While the pj spectra of the <f> in the whole pj range and of the 
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<jd at p T > 0.8 GeV/c are perfectly thermal and agree very well with the "blast-wave" fits at all 
centralities, the low-p T a> spectra become gradually depleted as centrality increases, with almost 
complete disappearence of the a> with pj < 0.2 GeV/c for the most central collisions. 

Significant attention was paid recently to so called 0-puzzle: the disagreement between 
4> -> K + K~ measured by NA49 E2 and (p -> ^ pT measured by NA50 03]. NA50, whose 
<p acceptance is limited by pj > 1.1 GeV/c sees nearly twice the yield observed by NA49 whose 
statistics is limited by p T < 1.6 GeV/c, with a significant difference in the inverse slope pa- 
rameters: 7^ = 234 + 7 vs. Tkk = 305 + 15 MeV in the central Pb-Pb collisions. Recent 
measurements by CERES |44| of (f> production in central Pb-Au collisions at 158 A GeV both 
in the K + K~ and e + e" channels are consistent with each other and seem to confirm the NA49 
results (although with large errors on the T e $ in e + e~). NA60 has reported similar measurements 
of both K + K~ and p + pT channels in In-In collisions [45 1. Like CERES, it finds a good agreement 
between the two channels, with T^iFig.^right)) compatible with the observations of NA49 for 
the same number of participants, but it observes a smaller <p (Fig. ^/e/fj) yield per participant 
than NA50 and slighly stronger than both NA49 and CERES do. Taking into account that NA50 
recently reanalysed its Pb-Pb data and confirmed the previous results [46], at present, it is dif- 
ficult to reconcile all of the observations into a coherent picture, albeit there is some hint for a 
possible physics mechanism leading to a difference in the two channels H71 . 
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3. Conclusions 

A summary of dilepton measurements at the SPS was presented. The most prominent result 
is the excess observed in heavy ion collisions at all masses below the J /if/. Its most plausible 
explanation by the production of thermal dimuons is supported by the thermal-like spectra both in 
mass and in transverse momentum, the lack of any polarisation and by the reasonable agreement 
with theory. At low masses, the l + T pairs are dominated by nn annihilation via the short-lived 
p meson broadened by the hot medium (the Brown-Rho scenario of a dropping p mass is ruled 
out by the data). The monotonic rise of r e ff with mass up to M ~ 1 GeV/c 2 with a sudden drop 
and then stabilization at higher masses, seen by NA60, suggests that in the IMR the dominant 
contribution is due to the qq anniliation in the partonic phase, when the radial flow has not yet 
developed. The suppression of the low-/?r ai mesons in the central In-In collisions hints on the 
first observation of its in-medium modifications. The 0-puzzle: the contradiction between the 
and K + K~ decay channels studied by the NA50 and NA49 experiments in Pb-Pb collisions 
is not solved despite the consistent results obtained for both channels by CERES and NA60 in 
Pb-Au and In-In collisions, respectively. 
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